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1. Introduction 

The phase transition in Nf = 2 QCD (m„ = m^/ = m) can provide fundamental insight into the 
mechanism of confinement[|l|]. 

In the limit m — > oo the system is quenched, the deconfining transition is first order, (L) is the 
order parameter and Z3 the symmetry. However the coupling to quarks breaks Z3. 

At m ~ the deconfining transition seems to coincide with the chiral transition. A transition 
line exists at intermediate values of m which is defined by the maxima of a number of susceptibili- 
ties (the specific heat Cy, the susceptibility Xm of the chiral condensate {^^1/), the susceptibility of 
(L)) which all coincide whithin eiTors (see e.g. Ref.[Q]). 

At m ~ a renormalization group analysis [|3|] indicates that the chiral transition is first order 
for Nf > 3, and for Nf = 2 can be either first order or second order in the universality class of 
0(4). For the latter case the transition is a crossover at m 7^ 0, and in particular- a tried tical point is 
expected in the T — [i plane [Q] which could be observed in heavy ion collisions. 

Various groups have studied this problem with Wilson[||] or staggered[^, ^ ^ |8|, ^ fermions. 
No clear sign of discontinuities is found at least for the lattice sizes used, but no agreement either 
with the critical exponents of 0(4). 

We investigate the issue using standard Kogut-Susskind fermions on lattices 4 x with Ls = 
12, 16,20,24,32 and lattice quark masses aniq ranging from 0.01335 to 0.307036. In the present 
work we give our final results using our whole dataset collected on APEmille machines and a 
detailed analysis which is different in many points from the previous ones present in the literature. 
Preliminary reports of this work were presented at past conferences. 



2. Analysis of the critical behavior 

The order of a phase transition can be investigated by finite size scaling analysis. Near the 
critical point, for a second or a weak first order transition, the singular behavior is described by 
power law divergences according to universal critical exponents. The singular part of the free 
energy density is a homogeneuos function of the reduced temperature z = 1 — T /Tc, the lattice 
quark mass aniq and the size of the system L^: 

^,,{T,amg,L,) = L;"^^, {xU; ,am^l7,") (2.1) 

There are two independent critical exponents: the thermal jt and the magnetic one yu. The other 
more familiar- critical exponents are given by (<i = 3 is the dimensionality of the system): 

a = 2-^- p = ^-'JL 7 = 2^-^ 5 = --^ v = i (2.2) 

yt yt yt yt yt d-yn y, 

The numerical values of the critical exponents of interest for this work are given in Table |I[ 



According to Eq. |2JJ this problem involves two scaling variables. We simplify the problem 



restricting the parameter space in three different ways: 

1. keep aniqLs'' fixed while varying L^. One must choose a fixed value for j/, so that only one 
universality class at a time can be tested. Note that by chance the numerical value of j/, is 
the same for 0(4) and 0(2), so in practice we can test these two critical behaviors at once. 
The scaling law for ^5 is given by: = L7^^, (tL?) where now aniq is a function of L^; 
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2. keep T and aniq fixed and take the limit oo (finite mass scaling, Ref. |8[, This 
amounts to assume the existence of one single scale {aniq) in the problem, and no divergence 
of the coiTclation length at the transition: i.e. finite mass scaling and no finite size scaling 

: Cirri, 



Eq. (gjj) becomes: ■ 



^si'^aniq^'^^'') where the size of the system is not present; 



3. the last possibility is to keep tIZ' fixed while taking much bigger than the pion cor- 
relation lenght: aL^ » ^/mjt- This should work better than case 2. if the con^elation 
lenght becomes large and comparable to L^. The free energy scaling law in this case is: 

In cases 2. and 3. we are not forced to fix j/, so we are free to test all the critical behaviors of 
interest. The price for this is that we were forced to make additional approximations. 

Taking appropriate derivatives we find the relevant scaling laws for the thermodynamic sus- 
ceptibilities^ . In this work we concentrate on two of them: the specific heat Cy and the chiral 
susceptibility Note that while the first quantity is always relevant, the second one is a good 
quantity only if (v/^i//^) is the order parameter: that is certainly true for m = but may not be for 
m > 0. The expected behavior for the cases listed above are: 



1. Cy — Co = l!^'^ ^d'^Ll'^) ; Xm — Xo = lJ^^ ^^{zLy^). i.e. the peaks scale as a power of L 



l/v. 



and the widths as L., 



l/v. 



2. Cv-Co = (am,)-«/^>'"0,(T(am,)->'/>''O ; Xm - Xo = (am,)-r/^>''.,^;^(T(am,)->''/>'"). i.e. the 
peaks scale as a power of amq and the widths as arriq 



3. Cy - Co = (am^)-«/^>'"0,(TLy'') ; Xm - Xo = (am^)-T'/^>'"0;^(TLy^). i.e. the peaks scale 



as 



a power of antq and the width as Lg 



-l/v 



To test different scaling hypotheses we have performed a huge amount of Monte Carlo sim- 
ulations, divided into three groups (see Table 0). The first two are called Run! and Run2 and 





yt 


yh 
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0(4) 


1.336(25) 


2.487(3) 


0.748(14) 


-0.24(6) 


1.479(94) 


0.3837(69) 


4.852(24) 


0(2) 


1.496(20) 


2.485(3) 


0.668(9) 


-0.005(7) 


1.317(38) 


0.3442(20) 


4.826(12) 


MF 


3/2 
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2/3 
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1/2 
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V Order 


3 
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1/3 


1 


1 
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Table 1: Critical exponents of (9(4), (9(2), mean field and first order. 





Runl 


Run2 




12 


16 


20 


32 


12 


16 


20 


32 


aniq 


0.153518 


0.075 


0.04303 


0.01335 


0.307036 


0.15 


0.08606 


0.0267 


#Traj. 


22500 


87700 


14520 


14500 


25000 


131390 


16100 


15100 


aLs ■ 


11.9 


11.0 


10.0 


8.9 


11.3 


15.8 


14.8 


12.4 



Table 2: Run parameters for our Monte Carlo simulations. A more detailed list can be found in Ref. [|l|] 



See Ref. nih for the detailed definition of the susceptibilities on the lattice. 
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Figure 1: Pseudocritical coupling values as a function of aniq. The fit (explained in the text) has /dof - 
0.9, dof = 15. Values obtained from other collaborations are also included. 



they satisfy the requirements of case 1., namely we have fixed aniqLl'' to two different values and 
yh appropriate for 0(4) and 0(2). The third group is formed by two other simulations with the 
following parameters: Ls = 16, amq = 0.01335; = 24, aiUq = 0.04444; with 10000 trajectories 
each. Simulations were made using the hybrid R algorithm sweeping a wide range of values of j8 
for each value of aniq. The datasets collected were analysed using the multi-histogram reweighting 
procedure. 



3. Results 

On the lattice the reduced temperature is expressed as a function of the parameters j8, amq. 

Tq a[p,amq) 

where a{^,amq) is the lattice spacing. Expanding a in power series near (j3o,0) we obtain: 

Toe {^Q-^)+k,„amq + k„^i{amqf + k„,pamq{^o-^) + --- (3.2) 

In previous works in the literature it was assumed that T (jSq — j8 ) thus neglecting the lattice 
spacing dependence on the lattice quark mass aniq. The terms in Eq. ( |3.2| ) prove to be sufficient 
to describe the data of our simulations. There are two different way the pseudocritical temperature 
could scale: = k^L\^^ or = k'^antq^''' . The expansion parameters and or k'^ can be obtained 
by fitting numerical data. The results are the following: there is no visible shift of Tc varying Lg, 
i.e. k^ is always compatible with zero; data are equally well described by all the critical behaviors 
expected (first order, 0(4), 0(2), MF) if one includes in the expression of Tc the quadratic terms 
shown above; including only the linear terms the statement doesn't change if we restrict the allowed 
range of atUq (see Ref. [jl]]); taking oc j8o — j3c does not describe the data in any mass range at our 
disposal for a second order while for a first order it is identical to keep only the linear mass term in 
Eq. (|3). 

Having fixed the parameters for the reduced temperature T (see Fig. |l]) we can next consider 
the scaling of Cy and Xm- We first consider case 1: using data from Runl and Run2 we can 
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Figure 2: Scaling behavior of Cy and Xm as expected in case 1. No sign of scaling is found. This excludes 
a second order (9(4) or (9(2) transition. Simulations from Runl are shown; Run2 is similar[|lj. 

check 0(4) and 0(2) critical behaviors. If there is scaling all the curves obtained by plotting the 
susceptibilities divided by the appropriate power of as functions of the scaling variable tL? 
should collapse onto a universal scaling function. As shown in Fig. ^ no such scaling is observed 
for 0(4) (0(2) is quite similar). This clearly excludes these universality classes. To check other 
critical behaviors we consider the scaling of cases 2. and 3. We show in Fig. ^ the analysis of Cy '- 
similar considerations also apply to ;^„,. The second order critical behavior does not describe our 
data. The growth of the peaks is instead consistent with a first order. This is an indication for a 
first order transition. The widths of the curves seem to be described by the scaling of case 3 and 
not according to case 2. This means that the dependence on cannot be neglected. We have also 
explicitely observed consistent deviations from the scaUng predicted in case 2. making simulations 




Figure 3: Scaling of Cy according to case 2. (top) and 3. (bottom) for both (9(4) (left) and first order (right). 
Other second order behaviors are similar to the one shown. 
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Figure 4: Scaling of (v^^V'^) for (9(4) (left) and first order (right). 



at two different volumes with the same aniq, see Ref. |[T]]. In fact this is due to the presence of two 
different physical scales in the problem. 

One can also study the magnetic equation of state, i.e. the scaling of {^^f) — {Y¥)o = 
amq^ F {x I aniq^^'') as shown in Fig. ^. Also in this case the first order behavior seems to de- 
scribe well the data while the second order is excluded. 

We have looked for metastabilities in the time histories as required by a first order transition: 
no clear evidence exists up to the lattice sizes explored. 

4. Conclusions and outlook 

We have investigated the nature of the chiral transition mNf = 2 QCD. By using the correct 
definition of T on the lattice we have shown that it is not possibile to discriminate the order of the 
transition by looking only at j8f with the present data. Studying the critical behavior at fixed aniqLl'' 
we were able to exclude the second order universality classes, i.e. 0(4), 0(2), mean field. We 
found evidence for a first order transition looking at the scaling of thermodynamic susceptibilities 
but no clear sign of discontinuities. The magnetic equation of state also seems compatible with 
a first order transition. We plan to repeat the analysis at fixed aniqL^'' for a first order transition 
employing an improved action and = 6 to reduce possible artifacts. 
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